Quantum localization and bound state formation in Bose-Einstein condensates 
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We discuss the possibility of exponential quantum localization in systems of ultracold bosonic 
atoms with repulsive interactions in open optical lattices without disorder. We show that expo- 
nential localization occurs in the maximally excited state of the lowest energy band. We establish 
the conditions under which the presence of the upper energy bands can be neglected, determine 
the successive stages and the quantum phase boundaries at which localization occurs, and discuss 
schemes to detect it experimentally by visibility measurements. The discussed mechanism is a par- 
ticular type of quantum localization that is intuitively understood in terms of the interplay between 
nonlinearity and a bounded energy spectrum. 
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1. INTRODUCTION 

The phenomenon of Anderson localization in disor- 
dered quantum systems [I| was originally discovered in 
the context of the study of electrons in a crystal with 
imperfections ,2j. In fact, it is much more general |3| and 
has been observed in a variety of systems, including light 
waves in random media [J,|5|. Despite remarkable efforts, 
Anderson localization has not been observed directly in 
crystals, owing to the high electron-electron and electron- 
phonon interactions. It has finally been observed in non- 
interacting Bose-Einstein condensates in one-dimensional 
quasi-periodic optical lattices [6] , that feature a crossover 
between extended and exponentially localized states, as 
in the case of purely random disorder in higher dimen- 
sions; moreover, the effects of random disorder in optical 
lattices can also be simulated manipulating the interac- 
tions in multi-species mixtures [TJ. These achievements 
are due to the unprecedented degree of control over the 
system physical parameters, in particular the vanishing 
of the interaction strength, that ultracold atoms offer. 

Indeed, ultracold degenerate gases in optical lattices 
provide an unprecedented toolbox for the experimental 
realization of what were once just toy models sketch- 
ing the key features of complex condensed matter sys- 
tems. One prominent example is the Bose-Hubbard 
model [1, [0] , that was originally introduced as a variant 
of the better known Hubbard model and whose prop- 
erties were later discussed at length in the context of 
the description of superfluid "^He trapped in porous me- 
dia. The suggested realization in optical lattices loaded 
with ultracold bosonic atoms [1^ was soon achieved in 
a spectacular breakthrough experiment 'lli|. Driven by 
this brilliant result, a growing number of investigations 
has focused on the possibility to use optical lattices to re- 
alize various phenomena of considerable interest in con- 
densed matter physics [l^, [Ijl- Amongst these, in the 
last years much attention has been devoted to the study 
of localized quantum phases in many-body systems. For 



instance, it has been showed that it is possible to use 
boundary dissipation [l4'| or the control of the sign of 
the local interactions, exploiting Feshbach resonances, to 
switch from the repulsive Hubbard model to the attrac- 
tive one, whose ground state may feature a collapse of all 
the atoms of the system into a single site of the lattice 
|15l - [l7| . The transition to collapse is essentially due to 
the combination of the nonlinear dependence of the local 
Hamiltonian on the site occupation that makes energet- 
ically favorable those states that are characterized by a 
concentration of all atoms in a single site. 

In the present work we describe a route to quantum 
localization in many-boson systems with repulsive inter- 
actions, that has one important feature in common with 
the transition to collapse discussed in Refs. |15l - [l7| . The 
basic idea is to consider the maximally excited state in 
the lowest energy band of an interacting system in a 
non-translationally invariant lattice. Within the Bose- 
Hubbard framework this is just the eigenstate with high- 
est energy of an Hamiltonian H with repulsive interac- 
tions. This state is then the eigenstate with lowest energy 
(i.e. the ground state) of a new Hamiltonian H equal 
to minus the original Hamiltonian: H = — iJ. The ro- 
tated Hamiltonian has attractive interactions (instead of 
repulsive interactions) and a negative tunneling ampli- 
tude. However, this can be turned positive again by a 
TT-phase shift on every other lattice site. If one can show 
the occurrence of exponential localization in the high- 
est excited state of the repulsive Hamiltonian, this phe- 
nomenon is then completely equivalent to the collapse in 
the ground state of the corresponding attractive model. 
In the latter case, upon increasing the intensity of the at- 
traction, the particles form a bound state, with increas- 
ing mass, which appears to be localized if the correlation 
length becomes smaller than the lattice spacing. The two 
mechanisms differ in that the former is not realized as an 
on-site collapse in the ground state of a system with at- 
tractive interactions, but rather as a proper exponential 
localization in the maximally excited state of the lowest 



energy band in systems with repulsive interactions. It 
is thus a mechanism that is solely due to the interplay 
between nonlinearity and an energy spectrum bounded 
from above. The latter in turn is a fundamental feature 
associated to the fact that the Bosc-Hubbard Hamilto- 
nian preserves the total particle number. 

In the following we will investigate the properties of 
the maximally excited state of the one-dimensional repul- 
sive Bose-Hubbard model defined on open lattice chains. 
Such analysis, besides resorting to exact diagonalization 
for systems of very small size, will be carried out for larger 
systems using numerical solutions obtained with a con- 
trollable recursive algorithm as well as a semi-classical 
approach when the size of the problem makes the nu- 
merical computation impractical. We will show that, de- 
pending on the physical parameters of the system, the 
maximally excited state of the repulsive Bose-Hubbard 
Hamiltonian defined on an open chain features three dif- 
ferent phases of which the first one, associated to small 
values of the local repulsion, is characterized by a rela- 
tive atomic population spread on all sites of the lattice. 
At intermediate values of the on-site interaction there 
occurs a second phase in which a macroscopic fraction 
of the atoms begins to localize in a single site while the 
remainder of the atomic population is still spread over 
the lattice. Finally, we will show that for large values of 
the on-site repulsion the maximally excited state is char- 
acterized by an exponential localization in the center of 
the lattice and we will investigate the decay rate both nu- 
merically and analytically. We will then determine the 
physical conditions such that the overlap of the maxi- 
mally excited state of the lowest energy band with the 
lowest states of the upper bands can be neglected, and 
we will discuss how to detect experimentally the three 
different behaviors by measurements of the visibility. 



2. MODEL AND METHODS 

Let us consider a system of N ultracold atoms with re- 
pulsive on-site interactions described by a Bose-Hubbard 
model on a one-dimensional lattice of M sites: 
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One needs to consider open chains to look, even in prin- 
ciple, for the possibility of true localization. Indeed, in 
a translationally-invariant geometry the atoms would be 
unable to localize on a definite site. Namely, even in the 
presence of strong repulsive on-site interactions, the max- 
imally excited state would be essentially a Schrodinger- 
cat state, i.e. a superposition of localized states charac- 
terized by a flat distribution of the atomic density over 
the entire lattice [l7J. In Eq.© d = {M - l)/2, a^ [a]) 
are the bosonic annihilation (creation) operators on the 
j-th site, fij = d^ttj are the occupation number opera- 
tors, [/ > is the strength of the repulsive nonlinear 



on-site interaction, and T is the hopping amplitude be- 
tween neighboring sites. 

In order to determine an optimized analytic approx- 
imation to the maximally excited state of the Hamil- 
tonian Eq.([l]) on a finite open chain we proceed by a 
dynamical variational method and compare results with 
the ones obtained by exact diagonalization. We follow 
the route adopted for the corresponding attractive model 
[la . Il7| , introducing a macroscopic trial state of the form 
\(t>) — e^'^lcf)) where ip is a time-dependent phase and \(j>) 
is a coherent state of the form 
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Here jfi) is the vacuum state and the coherent-state con- 



stants, 
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, M, must satisfy the normal- 
ization condition X^i^i \4>j\'^ — 1- The complex quanti- 
ties 0j describe the on-site bosonic states by the on-site 
population \4>j\^ and the macroscopic local phase aigcjij. 
The request that the trial state satisfies the Schrodinger 
equation on the average, {(l)\idt — H\(j)) = 0, identifies the 
time derivative </? with an effective Lagrangian for the 
dynamical variables 4>j and the corresponding effective 
Hamiltonian [isll: 
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Maximizing the latter with respect to the variables 4)j 
under the normalization constraint, one obtains a semi- 
classical variational approximation to the maximally ex- 
cited energy eigenstate of the system. The first term of 
Eq.Q, i.e. the on-site interaction term, does not depend 
on the phases of the on-site variables (f)j , while the hop- 
ping term, at arbitrarily fixed values of \(j)j \ and \(l)j+i \ , is 
maximized by a phase difference ±7r. Therefore, except 
for an irrelevant global phase factor, the values (jjj associ- 
ated to the maximally excited state can be assumed to be 
real quantities with alternating signs. Defining Xj = \(t>j\, 
and taking into account the property of invariance under 
mirror reflection (xj — x^j) with respect to the center 
of the flnite chain, Eq.Q can be recast in the equivalent 
forms 

^ = |(A^-1)(4 + 2E4|+2TE^.-^.-i. 




Introducing the dimensionless ratio of the interaction to 
kinetic energy scales: A = U(N — l)/(8r), one can solve 
the problem for different values of A and maximize Eq. @ 
using the hyperspherical representation of the variables 
Xj to enforce automatically their exact normalization. 



A=0.05 




FIG. 1: (Color online) The atomic density profile {nj}/N in 
the maximally excited state of Hamiltonian ((T} obtained by 
exact diagonalization (full black squares) and by the semi- 
classical variational approximation (empty red circles linked 
by solid red line), for different values of the dimensionless 
energy ratio A in an open chain of M = 9 sites and A*' = 9 
atoms. All quantities being plotted are dimensionless. 



and can then compares the atomic distribution densi- 
ties so obtained with the ones provided by exact diag- 
onalization for small samples. The latter in turn can 
be performed very efficiently with the help of augmented 
recursive Lanczos algorithms [l^ . The result of this com- 
parison is reported in Fig.(IT]). 

Fig.© shows that the semi-classical solution provides 
an excellent approximation to the maximally excited 
state that becomes more and more accurate with in- 
creasing strength of the interaction and of the localiza- 
tion of the atoms at the center of the lattice, and can 
thus be extended to systems of much larger size that 
cannot be investigated by exact diagonalization. More 
important, Fig.([T]) shows that as the A parameter varies 
the atomic density profile in the maximally excited state 
crosses three different phases. The first one, associated 
to small values of A, is characterized by the absence of 
localization; the third one, associated to very large val- 
ues of A, corresponds to a complete concentration of all 
the atoms of the system in the center of the lattice; fi- 
nally, the second one, associated to intermediate values 
of A, corresponds to the onset of the localization of a sig- 
nificant fraction of the atomic population in the central 
site, while the distribution of the remainder of the atomic 
population over the entire lattice stays finite. 




FIG. 2: (Color online) Atomic density distribution (in log- 
arithmic scale) associated to the maximally excited energy 
eigenstate in the localized phase obtained by exact diagonal- 
ization (solid lines and full symbols) and by the semi-classical 
variational approximation (dashed lines and empty symbols) 
for a chain of Af = 9 sites and A'' = 9 atoms. From top to bot- 
tom: Black circles and lines (A = 1.0); red squares and lines 
(A — 2.0); blue diamonds and lines (A — 3.0). All quantities 
being plotted are dimensionless. 



EXPONENTIAL LOCALIZATION AND 
VISIBILITY 



We now investigate in detail the exponential nature 
of the localization in the maximally excited state. The 
presence of an exponential decay is reported (in loga- 
rithmic scale) in Fig. ([2]). It shows the behavior of the 
relative occupation {rii) /N as a function of the distance 
from the center of the lattice. The occurrence of an ex- 
ponential localization allows to introduce a simple and 
effective method to obtain an excellent analytical ap- 
proximate solution the problem of the maximization of 
the effective Hamiltonian. We introduce a dimension- 
less parameter x ^nd assume |0j| = x' ^j G [^id], with 
|0o| = 1 ~ 2x^(1 — X^'')/(l — X^) in order to satisfy nor- 
malization. Within this setting the maximization of the 
effective variational Hamiltonian can be performed ana- 
lytically. Exploiting the condition x ^ 1, necessary to 
have a state localized in the central site, one finds the 
approximate analytical solution x = A/8. It is worth 
observing that the expression for x that determines the 
approximate solution to the maximally excited state of 
the system does not depend on the size of the chain. On 
the other hand, it depends on the number of atoms in 
the lattice through A. In Fig. ([3]) we have compared the 
exact numerical solution of the variational problem with 
the approximate analytical solution for different value of 
A within the interval compatible with localization. One 
sees that even at moderate values of A the analytical ap- 
proximation reproduces the essential features of the exact 
localized quantum solution. 

The onset of localization in the maximally excited state 
can be naturally captured either by analyzing the behav- 
ior of the relative occupation of the central site {no) /N 
or by looking at the factor /o = 1 — {riQ) /N that mea- 
sures the relative atomic population outside the central 
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FIG. 3: (Color online) Atomic density distribution (in loga- 
rithmic scale) as a function of the distance from the center 
of the lattice for an open chain of M = 21 sites and A^ = 21 
atoms. Solid lines and filled symbols: Exact numerical so- 
lution of the variational problem. Dashed lines and empty 
symbols: Approximate analytical solution x = A/8. From 
top to bottom: Blue lines and triangles (A — 2.0); black lines 
and circles (A — 5.0); red lines and squares (A — 10.0); green 
lines and diamonds (A = 20.0). All quantities being plotted 
are dimensionless. 
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FIG. 5: (Color online) Relative atomic population outside 
the central site /o (red squares) and relative population of 
the upper energy band r (black circles) as a function of the 
gap width e for different values of the energy ratio A. Filled 
symbols: A = 2. Empty symbols: A = 3. All quantities being 
plotted are dimensionless. 



the momentum interference patter n l6|| . The visibility 
is related to the overall coherence [20| according to the 
relation 
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where Smax and Smm are the maximum and minimum 
values of the momentum distribution function: 
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In the lower panel of Fig.((4l) we report the behavior of 
the visibility, as a function of A. We see that the visi- 
bility is very sensitive to the onset of localization. It is 
approximately constant around its maximum when the 
atoms are delocalized over the lattice at small values of 
A and begins to decrease exponentially, at the onset of 
the transition in the center of the lattice, for Ac ~ 0.7. 



FIG. 4: (Color online) Relative atomic population ratio /o — 
1 — {no)/N (upper panel) and visibility V (lower panel) as 
functions of A for open chains of different size at unit filling. 
Green solid lines: Chain of M — 11 sites and A'^ = 11 atoms. 
Red dashed lines: Chain of M — 21 sites and A = 21 atoms. 
Black dot-dashed lines: Chain of M — 31 sites and A = 31 
atoms. All quantities being plotted are dimensionless. 



site. The upper panel of Fig.(j4l) shows that for small 
values of A the ratio /o is enhanced. In this regime and 
for very long chains, such that the border effects can be 
neglected, /o — >■ 1 — 1/M. It begins to decrease at the 
onset of localization at the critical value A^ ~ 0.7, finally 
vanishing asymptotically in the limit A — ^ 00. 

The rationale for the study of /o lies in the fact that 
it is directly associated to the visibility V, a quantity 
that can be actually measured by looking at the relative 
difference between the maximum and the minimum of 



4. ROLE OF HIGHER BANDS 

A most serious issue concerns the possibility that in 
the process of driving the system in the maximally ex- 
cited state of the lowest energy band, the lowest states of 
the higher energy bands may get significantly populated. 
The problem is then to establish under what conditions 
the relative overlap between these states and the maxi- 
mally excited level of the lowest band is negligible. We 
thus need to consider the two-band Bose-Hubbard Hamil- 
tonian with intra-band and inter-band interaction terms 
(21| . One needs to add to the lowest-band Bosc-Hubbard 
Hamiltonian Eq. ([Ij both the Bose-Hubbard Hamiltonian 
of the first upper energy band 
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and the inter-band interaction terms 
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In Eqs. (|7l8p Eg is the energy gap between the first ex- 
cited level of the optical lattice potential and the relative 
ground state, while d^- , dj- , and n^ are, respectively, 
the on-site bosonic annihilation, creation, and number 
operators relative to the first upper energy band. The 
Hamiltonian parameters of the two bands are obviously 
not independent: resorting to the standard harmonic ap- 
proximation for the optical lattice potential one has that 
T2 ~ 9.4T, C/2 ~ 3C//4, and W ~ U/2. Hence, the to- 
tal Hamiltonian Ht = H + H2 + Hj depends only on 
the two independent parameters e = Eg/8T, that mea- 
sures the width of the energy gap, and the previously 
introduced A = U{N — 1)/(8T), that expresses the ratio 
of the interaction to kinetic energy scales. Going again 
through the same dynamical variational procedure for 
the two-band model Ht and finite values of the gap e, 
we solve the maximization problem in the previously de- 
termined range of values of A that are compatible with 
exponential localization in the maximally excited state 
of the lowest energy band. As we can see from Fig.®, 
the relative atomic population outside the central site /o 
is quite unaffected by the inter-band energy gap width 
until the latter becomes so small to allow a relative pop- 
ulation r — J2j=-di''^'j)/^ of the first upper band that is 
comparable with /q. If the energy gap is further reduced, 
the occupation outside the central site begins to increase 
exponentially even if the population of the central site re- 
mains a substantial fraction of the total number of atoms 
until £ ~ A and Eg ~ U{N — 1). This finding allows to 
conclude that if the lattice is loaded with a total number 
of atoms N < Eg/{10U), one can safely disregard the 
presence of the upper energy bands. 



5. CONCLUSIONS AND OUTLOOK 

In conclusion, we have introduced and discussed a 
mechanism of exponential localization in the maximally 



excited state for systems of ultracold bosonic atoms with 
repulsive interactions in open optical lattices. The prop- 
erties of the maximally excited state have been studied as 
a function of the Hamiltonian parameter both with nu- 
merical and analytical techniques in order to determine 
the region of the parameter space in which exponential 
localization take place and the dependence of the expo- 
nential decay on the Hamiltonian parameters. Finally, 
we have discussed how the transition to localization can 
be detected experimentally by visibility measurements, 
and we have established the physical conditions under 
which the overlap with the upper energy bands can be 
neglected. This localization mechanism depends on the 
properties of the maximally excited state that stem from 
the interplay of nonlinearity and a bounded energy spec- 
trum. Being not a ground state property, it does not 
require the presence of random disorder, attractive local 
potentials, or ad hoc truncations of the Hilbert space. 

At first sight it would seem that to populate the max- 
imally excited state of a system in realistic experiments 
is an extremely challenging goal to achieve. However, 
thanks to the exceptional properties of controllability 
and manipulability of optical lattice systems, the pre- 
dicted phenomenon might be observed first by cooling 
the system in the presence of a strong local field favor- 
ing a substantial atomic population at the center of the 
lattice. After switching off instantaneously the local field 
the system would remain, with probability close to unity, 
in the strongly localized maximally excited state. Re- 
peated transitions to delocalization and exponential re- 
localization could then be observed simply varying the 
energy ratio A by changing the depth of the lattice and/or 
by tuning the scattering length. 
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